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Abstract 

The effects of space charges on hysteresis loops and field distributions in ferroelectrics have been 
investigated numerically using the phenomenological Landau-Ginzburg-Devonshire theory. Cases 
with the ferroelectric fully and partially depleted have been considered. In general, increasing the 
number of charged impurities results in a lowering of the polarization and coercive field values. 
Squarer loops were observed in the partially depleted cases and a method was proposed to identify 
fully depleted samples experimentally from dielectric and polarization measurements alone. Un- 
usual field distributions found for higher dopant concentrations have some interesting implications 
for leakage mechanisms and limit the range of validity of usual semiconductor equations for carrier 
transport. 
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I. INTRODUCTION 



While in bulk they behave like insulators, the description of ferro electrics as wide-band- 
gap semiconductors is now widely used to account for many properties of thin films. The 
semiconducting state is, in general, thought to be detrimental to ferroelectric properties and 
thus adaptation of standard semiconductor theory to ferroelectric materials has been the 
theme of many recent publications. However, despite numerous theoretical efforts, no con- 
sensus has yet been reached on how to include the effects of semiconductivity when analysing 
experimental data. Analytical models tend to make simplifying assumptions which, while 
providing intuitive insight, may miss out important effects or restrict the range of validity 
of the model to less interesting cases from the experimental point of view, whereas detailed 
numerical simulations often assume many parameters that cannot be easily extracted from 
experiments. A common analytical approach is to include a uniform ferroelectric polariza- 



tion as two sheets of opposite charge diplaced from the electrodes 



dv some distance assumed 



to be the thickness of a non-ferroelectric (or "dead") layer |2j, y|. Pintilie et al. 
have used such a model to interpret their electrical measurements on semiconducting epi- 
taxial single crystal (Pb,Zr)Ti03 (PZT). However, as we will show later, the assumption of 
a uniform polarization may become invalid when space charge levels exceed ~ 10 19 cm~ 3 , 
and spatial variation of polarization P(z) has to be taken into consideration. 

One way of including a non-homogeneous polarization is in the framework of Landau- 
Ginzburg-Devonshire (LGD) theory. The free energy density is expanded in powers of P{z) 
and the gradients VP(r) truncated at an appropriate power; its minimization using the 
variational principle generally results in a second order differential equation whose solution 
requires the boundary conditions to be known. Since no agreement has been reached on 
what these boundary conditions are, they currently provide a playground for theoretical 
investigations into various surface effects, but at the same time they have a large influence on 
the polarization profiles within the bulk of the films, as was shown numerically by Baudry 0, 



Ignoring the gradient term in the free energy expansion leads to a polynomial equation, 
which was investigated analytically by Bratkovsky and Levanyuk They mainly looked 
at the effects of space charges on phase transitions showing a lowering of the transition 
temperature and suppression of ferroelectricity in the near-electrode layers. However, only 
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cases where the space charge was localised as either one central sheet or two symmetrically 
placed sheets of opposite charges were considered, and some of the analytic results are only 
valid for low doping concentrations. A polynomial equation was also addressed numerically 
by Brennan |8j], who included the effects of polycrystallinity and looked at uniform doping 
but again only for low doping concentrations. 

Finally, we would like to mention the work of Morozovska 111 Q who approached 

the problem of space-charge effects using the Landau-Khalatnikov equation, which has a 
time-dependent term and therefore allows the simulation of hysteresis loops at various fre- 
quencies. Apart from confirming the constriction of the hysteresis loops shown by Baudry 
and Bratkovsky and Levanyuk, Morozovska et al. have shown that dielectric losses may be 
responsible the tilting of the loops, as well as reported other effects such as the appearance 
of double hysteresis loops. Due to the complexity of the equations, however, the spatial 
variations of polarization and electric field were not investigated and only volume-averaged 
values were reported. 

None of the studies above (apart from 8]) have discussed in any detail the electric fields 
inside ferroelectric samples with space charge, and to date researchers of leakage currents 
in ferroelectrics are divided into those who interpret their data assuming a constant field, 
compatible with fully depleted semiconductors or insulators, and those who use the standard 
results for partially depleted semiconductors where the field drops linearly across a depletion 
layer to a negligible value in the bulk. 

We have focused on more realistic space charge distributions and simulated the electric 
field distributions and ferroelectric hysteresis loops in fully and partially depleted ferroelec- 
tric films using parameters comparable to those extracted experimentally for a high quality 
144nm thick PbZro.4Tio.6O3 thin film as an example. We show a decrease in coercive field 
and remanent polarization, P r , values with increased dopant levels as observed in previous 
works, but also some unusual field profiles and discuss their significance for leakage mech- 
anisms and coercive fields in ferroelectric thin films. Last, we discuss the main problem of 
the thermodynamic modelling approach, namely its overestimation of the coercive field Ec 
if bulk parameters are used. It will be discussed how the experimentally measured values of 
Eq can be reconciled with the phenomenological approach provided space charge is taken 
into account. 
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II. MODEL AND ASSUMPTIONS 



Landau-Ginzburg-Devonshire theory 

In LGD theory, the free energy F is composed of a bulk term F&, given as a power series 
expansion in P, and a surface term F$. The energy of a dipole p in a field E is — E ■ p; 
however, when considering the interaction of the local polariasation P(z) with the local field 
E(z) one has to be careful not to include the self-field, hence the factor of a half in front of 
the depolarization field (E^p) term in equation (where S is the cross-sectional area, E r 
is the electric field excluding the depolarization term and a, (3 and k are the usual Landau 
coefficients; A is the so-called extrapolation length Q|; a second order phase transition has 
been assumed). 

f = £dz \p\z) + \p\z) + I (^) 2 - E r (z)P(z) - \E dep {z)P{z) (1) 

^ = £ (P(0 + ) 2 + P(d-) 2 ) (2) 



F = F b + F s (3) 

Only terms up to fourth power in P(z) have been retained for the calculations in this 
paper, however, higher powers can be added without any additional difficulty when modelling 
first order phase transitions or bulk tetragonal PZT, which has high polarization values thus 
not allowing the P 6 (z) term to be neglected. 

Using the variational method to minimize the free energy yields a non-linear second order 
differential equation 

aP(z) + (3P 3 (z) - k^P- = E r (z) + P ~ P{Z) (4) 

CiZ 6g 



with boundary conditions determined by the surface term F t 



s 



fdP\ P(0+) fdP\ _Pidr)_ 

\dz) z=Q - X \dz) z=d X ■ 1 ) 



An approximate analytic solution to equation (J3J) was obtained by Glinchuk and Moro- 
zovska ^| for the case of an ideal insulating ferroelectric with no space charge. Baudry 



attempted to include the space charge effects by solving the equation numerically and found 



an interesting interplay between the surface effects due to the imposed boundary conditions 
and the electric field due to the space charge. However, due to the unspecified physical na- 
ture of the surface effect the calculated field and polarization profiles in his case seem to be 
inconsistent with the electrostatic equations used. The main problem in solving equation (J3J) 
is that the parameter A is unknown. Recent ab initio calculations for epitaxial KNbC>3 with 
SrRuOs and Pt electrodes by Duan et al. have shown that what happens near the 
surface is highly dependent on the type of atomic bonding between the ferroelectric and the 
metal, with both reduction and enhancement of local polarization possible. 

In this paper we are only interested in the effects of space charge and will therefore 
ignore what happens at the interface, especially since so little is known about the boundary 
conditions that should be used there. This assumption is most realistic for thicker films. We 
also neglect the gradient term in the bulk, which will be fully justified later. 

Electric fields 

The sample is assumed to be a uniformly doped monodomain single crystal ferroelectric 
of thickness d with symmetrical and perfectly metallic electrodes. The sum of the charges 
inside the sample and on the electrodes must be zero to have zero field in the electrodes. 
We have chosen to separate the charge distributions as illustrated in figure [TJ It can be 
shown that the polarization P(z) inside the film will induce a charge density a> on the 
electrodes equal to the average polarization P 



For a partially depleted sample o~ sc> x and <t SCj 2 will compensate the charges in the two de- 
pletion layers. In the case of full depletion at zero applied voltage a SCj i = a SCj i = 1/2 x total 
space charge per unit area in the sample; a ext and — a ext provide the external field. 
The electric field at any point in the sample is given by 




e E(z) = D{z)-P(z) 




o-ext + o- sc + qN D z + P - P(z) 




e E ext + D sc (z) + e E dep (z) 



(6) 
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FIG. 1: Volume charge densities (p sc and /?p) inside the ferroelectric and surface charge densities 
{(T ex t., a sc and op) on the electrodes assumed in the model. 

where D(z) is the electric displacement and E r (z) = E ext + D sc (z)/e . Inserting this into 
equation (j3J) and neglecting the gradient term leads to a third order polynomial 

aP(z) + [3P 3 (z) = E ext + + P ~ P{Z) (7) 

which, however, does not have a simple analytic solution due to the P term. This equation 
is the same as previously investigated by Bratkovsky and Levanyuk for the simple cases of 
a single cetral sheet of charge and two such sheets at symmetric planes. 



Strain effects 



In-plane strain caused by the mismatch between the lattice parameters of the ferroelectric 
and the substrate can either enhance or reduce the polarization of the ferroelectric compared 
to the bulk value depending on its sign. Previous works 12] show that the effect of 
strain is to simply rescale the local Landau coefficients a and (3, thus if we assume that the 
effective coefficients (i.e. those that would be extracted experimentally from the temperature 
dependence of the dielectric susceptibility and spontaneous polarization) are also simply 
rescaled, then our model should be valid for strained films too. 
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III. SIMULATION RESULTS 



A. Full depletion case 

In the case of full depletion, the fields entering equation (J7j) are given by 

p (Text _ Y_ 

-'-'ext — 7) 

eo a 
D sc (z) = - q N D {--z), 

assuming p sc {z) = qNp, as would be in the case of an n-type material; however, the model 
is equally valid for p-type samples with a simple sign reversal of the appropriate terms. 
Equation (J7J) was then discretized and solved numerically for various dopant concentrations 
Nd. The simulations were carried out as follows: an initial guess was made at the polar- 
ization profile, P was computed and equation (J7J) was solved to obtain the new polarization 
profile and thus a new P. The procedure was repeated until the convergence criterion 
(AP < 10 _11 C/m 2 for that iteration) was satisfied. The applied voltage was computed by 
integrating the total field in the sample. The following Landau coefficients were used in 
the simulation for the results to be comparable to our experimental data discussed later: 
a = -6.36 x 10 7 mF" 1 and (3 = 1.38 x 10 9 n^C^F -1 (giving a P s = y/-a/0 of 0.215 
Cm -2 ); d was assumed to be 150nm. 




Bias (V) 

FIG. 2: Effect of space charge on ferroelectric hysteresis loop for a fully depleted sample. A 
reduction from "ideal" to zero is observed in P r and E c over a very narrow range of Nd- 

Figure El shows that for concentrations below about 10 18 cm -3 the hysteresis loops remain 
unchanged from the ideal case of uniform polarization. At 10 19 cm -3 a clear shrinking of the 
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loop is observed and doubling this concentration leads to an almost total vanishing of P r . 
So within less than an order of magnitude in Nd the shape of the hysteresis loop changes 
from almost "ideal" to a simple curve. 
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FIG. 3: Spatial variation of (a) the local polarization P and (b) the total electric field E(z) = 
E sc {z) — (P — P(z))/eo inside the film at zero E ex t. Note, although P{z) is linear to a very good 
approximation, it is not exactly linear as can be seen from the curvature of the total field; this, of 
course, is no surprise, as equation (J7|) contains a cubic term. 

To understand these results let's look at the internal field distributions (figure EJ). The 
total field is given by 

E(z) =aP(z)+(3P 3 (z). (8) 

Due to the high dielectric constant of ferroelectrics, to a very good approximation D(z) = 
P(z). From Poisson's equation 



D{z) =P + D sc (z) + ^e . 



(9) 



The last term is generally much lower than either the first or the second for the usual 
applied fields and we will neglect it for simplicity. For low space charge concentrations 
the displacement is approximately constant throughout the film and therefore so is the 
polarization, perturbed only slightly by the linear D sc (z) i.e. P(z) = P + D sc (z) ~ P s + 
D 8C (z). Inserting this into (jHJ and neglecting terms of order D 2 sc (z) and higher we obtain 



E(z) 



(a + 3pP^)D sc {z) 
= -2a ■ D sc {z) 
= -2a ■ qN D (d/2 - z) 

i.e. a linear internal field profile just as in the case of a semiconductor with (ee ) _1 = —2a. 



S 



On the other hand, when the space charge dominates (f9*J). the polarization P(z) ~ D(z) ~ 
D sc (z) is approximately linear and antisymmetric about the midpoint of the film, averaging 
to zero. The total field given by (JSJ) is now 

E{z)^aD sc {z)+f3Dl c {z) (10) 

and is a cubic in z as can be seen in figure EJ This is very different from the usual field 
distribution is a semiconductor with space charge. 



B. Partial depletion 

For a semiconducting film with Nd = 2 x 10 19 cm~ 3 to be fully depleted at zero applied 
bias it must have a built in voltage Vn = qd 2 N£>/(8efeo) ~ 1.2V (assuming a uniform field 
independent dielectric constant e/eo = l/2a = 888eo), whereas the expected values are 
closer to a fewx0.1V thus for high dopant concentrations the films might be partially 

depleted. 

We will again assume a uniform density of fully ionised donor impurities. In the usual 
abrupt depletion approximation any region with a non-zero field is completely depleted of 
free charges and has positive space charge of volume density qNjj. In the regions free of 
space charge, the total field is zero 

E(z) = E ext --(P-P(z)) = 
e o 

and 



P(z) = P S = 




giving 

E ex t — — (P — Ps)- 

For simplicity we assume that all the voltage is dropped across the reverse biased depletion 
region of width w = w(V); then D sc (z) = ±qN[)(w — z) within the depletion region. This 
assumption is true for Vu = and is a good approximation for all but very low applied 
biases when Vu ^ 0. 

In the simulation, D sc (z) is set by choosing the depletion width w and the correspond- 
ing potential difference between the two electrodes is again calculated by integrating the 
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FIG. 4: Field profiles inside the 150nm film for different extents of depletion (i.e. at different 
biases) for four dopant concentrations: Np = 10 18 , 3 x 10 18 , 10 19 and 3 x 10 19 cm~ 3 . 

total internal field once the calculation converges and P{z) is found. The results are very 
interesting. 

We start with the sample poled positively (i.e. we force P{z) = +P S in the space charge 
free region) . Figure 0] shows that for small Njj the field increases approximately linearly 
with z as usual in the case of semiconductors. The ratio D sc (z)/eoE(z) tends to e with 
decreasing dopant concentration as expected. At larger Np, nonlinearities in E(z) due to 
the cubic term in equation ((7|) become obvious, causing an upward bending of the total field 
near the electrode. This has great implications for current- volt age characteristics. While 
the bias is increasing, the field near the interface is reduced and thus the usual field induced 
potential barrier lowering works in the opposite way - the barrier is raised, possibly leading 
to negative differential resistivity. This should be explored further. 

At still higher concentrations (figure EJ^c) and (d)) a point is reached when the field at 
the interface changes sign before full depletion is reached. The bias cannot be increased 
beyond this point without switching of the polarization. Changing P(z) to — P s in the 
undepleted regions gives us the rest of the negative half of the polarization loop. When 
Nd pa 3 x 10 19 cm -3 a second sign reversal is possible (figure EJl). 

The hysteresis loops for three doping concentrations are shown in figure 03 between —5 
and 5V. For the lowest concentration shown full depletion occurs at ±0.25V and the P-V 
curve shows a good overlap with the "ideal" loop. Increasing the doping leads to a reduction 
in the coercive voltage and the non-zero field polarization. However, P r remains the same 
as there is still no internal field at zero bias for Vy, = 0. In a way, the ferroelectric properties 
are improved with lower coercive fields, squarer loops and no reduction in the remanent 
polarization. 

For finite depletion regions form at both electrodes at zero bias, introducing space 
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FIG. 5: P-V loops for partially depleted samples between ±5V. The Np = 1 x 10 18 cm -3 loop is 
only shown in the range where the sample is partially depleted; outside this region the sample is 
fully depleted and is shown in figure HJ 




FIG. 6: (a) Variation of remanent polarization with built-in bias, (b) Built-in biases which can be 
supported without full depletion as a function of Nd for a ferroelectric and a semiconductor with 
eeo = —l/2a. Region X is not allowed for either, region SC is only allowed for the semiconductor, 
FE for ferroelectric and SC+FE for both. 

charge into the sample and thus lowering the remanent polarization in a similar way as in 
the fully depleted case. Figure a) shows how affects P r . As Nd increases higher voltage 
drops can be supported before full depletion is reached, as is usual in semiconductors, but 
within a certain range of Nd the electric fields in the deletion regions start to change sign and 
the maximum voltage drop that can be supported decreases. Part (b) of the figure compares 
the built-in voltages that can be supported by the space charge regions in a ferroelectric and 
a semiconductor of the same dielectric constant. The maximum around 10 25 m~ 3 corresponds 
to the point at which the electric field at the interface changes sign for the first time due to 
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the cubic term, whereas the minimum around 3 x 10 25 m -3 corresponds to a second reversal 
of sign. 

Stricly speaking, at very low biases, half of the applied voltage is dropped across each 
depletion region in a semiconductor. As the applied bias increases, the fraction dropped 
across the forward biased depletion region ff rapidly drops to zero. The exact fraction of 
the voltage dropped across each depletion width can be determined by equating the currents 
flowing through both regions, but this involves assuming a specific leakage mechanism. 
Usually, however, ff is only significant at very low biases, thus we can simulate hysteresis 
loops for samples with non-zero Vu assuming all the bias is applied across the reverse biased 
junction. In general, the higher the built-in bias the more space charge there is in the 
film and the more the loops will tend to those of fully depleted films with same dopant 
concentration. 



C. Comparison with experiments 

For an ideal bulk ferroelectric, whose free energy may be truncated at the P A (z) term, 
the a and j3 coeffcients can be found experimentally by measuring the zero field dielectric 
constant and remanent polarization P r and using the following relations 

eeo X e C 

where C is the Curie constant. Parasitic series capacintaces which reduce the effective 
measured dielectric constant can be eliminated by looking at the temperature dependence 
of e. Provided it is temperature independent, the parasitic contribution will not affect the 
slope of the 1/% vs. T plot which gives the Curie constant and will change only the intercept. 
Provided the effect of strain is to simply rescale the LGD coefficients then if these modified 
a' and (5' are extracted from zero-field measurements as discussed above, one expects that 



the hysteresis loop should fit the cubic equation [17 1 

a'P + (3'P 3 = E + «e (13) 
where £ is a constant and E = V/d. But what happens when we have space charge? 
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FIG. 7: Theoretical P — V plots for fully depleted films with dopant concentrations = 
10 18 ,10 19 and 1.5 x 10 19 cm~ 3 and the curves calculated using equation (|13|) with £ = and a' 
and (3' extracted using the simulated e(T) and P S (T). 

We have used our model with C = 4.44 x 10 5 K and T c = 548K to simulate the temperature 
dependence of P s and x, where \ was calculated from polarization values at the tickle voltage 
of 50mV (as it would be experimentally) using 



For partial depletion two cases were considered: 1) the applied 50mV drops entirely across 
the reverse biased contact i.e. ff = 0, and 2) half the applied bias drops across each contact 
i.e. ff = h the difference between the extracted Curie constants for the two cases was 
negligible. In what follows, a' and f3' are coefficients that would be extracted experimentally 
from e(T) and P r (T) measurements, modified from the bulk values a and (3 by strain and 
space charge. 

From figure [7| we can see that in fully depleted films the experimentally extracted a' 
and P', when substituted into equation (JTHJ), are expected to reproduce the hysteresis loops 
perfectly. This is not the case for partially depleted films. Figure Efa) shows that the 
experimentally extracted parameters do not give a good fit to the loop for all but very low 
fields; in fact the loops are not well described by the cubic (|T3|) even if we allow all the 
parameters to float in the fit. Decreasing the dopant concentration but keeping the built-in 
bias constant at 0.2V decreases the voltage range over which the film is partially depleted 
(figure Efb)). The hysteresis loop starts to resemble closer that of a fully depleted loop and 
the fit using experimentally extracted a' and f3' coefficients is expected to be better. 



X = 



P(50mV) - P(-50mV) 
lOOmV 



d. 
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FIG. 8: Theoretical P — V plots for partially depleted films with Vu = 0.2V and (a) Np = 
1.5 x 10 19 cm -3 and (b) Np = 10 19 cm~ 3 are shown as open cirles. The thick solid curves were 
calculated using equation 1)13(1 with £ = and a' and (5 1 extracted using the simulated e(T) and 
P S {T). Curves expected for fully depleted films with same are also shown for comparison 
(dashed). 

The extent of depletion in ferroelectrics has been a subject of debates for many years (for a 
discussion see chapter 5 of [18]), with some authors believing them to be almost always fully 
depleted due to their insulating nature, and others arguing them to be sometimes partially 
depleted due to the semiconducting behaviour often observed in thin films. Extracting 
the LGD coefficients using zero field dielectric and polarization measurements and looking 
at how well they recreate the observed hysteresis loops can in principle give a qualitative 
estimate of the extent of depletion in a film. In practice, other factors can affect the shape 
of the hysteresis loops. Parasitic series capacitances can cause a stretching on the loop along 
the field axis, increasing the apparent coercive fields. Morozovska et al. also showed that 
free carriers can result in tilting of the loops, which is likely to be significant in the case of 
semiconducting films. Thus, using the procedure described here, a good fit might indicate 
full depletion, whereas a poor fit may be caused by factors other than partial depletion. 

As an example, we present some experimental measurements on 144nm thick 
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PbZro.4Tio.6O3 capacitors from Samsung. The films deposted by a sol-gel method from 
99.9999% purity targets have a columnar structure on top of the [111] oriented Pt bottom 
electrodes, providing almost epitaxial interfaces for individual grains. Together with the 
Ir/Ir02 top electrodes the structures are fully integrated and embedded in Si. Further de- 
tails of the samples and their full electrical characterisarion will be published elsewhere, 
here we present only the temperature dependence of the capacitance measured with a 50mV 
tickle bias using the HP4192A impedance analyser (figure EJa)) arid the remanent polariza- 
tion obtained using the Radiant Precision Pro tester (figure Efb)). From the plot of e _1 (T) 
the Curie constant was calculated to be 4.44 x 10 5 K; P r 2 (T) plot yielded a value for f3' of 
1.52 x 10 9 m 5 C _2 F _1 and T c = 620K. A room temperature hysteresis loop is shown in fig- 
ure Etc) together with the LGD curve (fT3|) using the calculated parameters (and £ = 0). 
The excellent fit over the whole range of biases investigated indicates that the material 
is probably fully depleted [l^. This is consistent with the large resistivity expected and 
observed in these high purity samples, as well as with the leakage mechnism believed to 
be due bulk-limiting Poole-Frenkel conduction (the current- volt age data will be published 
separately) . 
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FIG. 9: Temperature dependences of (a) e , (b) P r 2 , and (c) LGD fit to the 1kHz hysteresis loop 
using a' and (5' extracted from (a) and (b). 

So far we have only discused how the extent of depletion affects the shape of hysteresis 
loops without mentioning the conditions necessary for full or partial depletion to occur. 
Figure ITU1 shows regions of full (shaded) and partial depletion as a function of donor impurity 
concentration and film thickness. For each region the expected shape of the hysteresis loop 
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FIG. 10: Regions of full (shaded) and partial depletion (unshaded) assuming Vbi = 0.1V, P r = 
20/iC/cm 2 and e = 890. 

is also sketched. Samples with a built-in voltage Vu and with a dielectric constant of e will 
be fully depleted at zero applied bias if 

8ee V bi 



N D < 



qd? 



(14) 



giving the line with a slope of —2 on the logiVo vs. log<i plot in figure Within the fully 
depleted region the P — E loops should be close to the "ideal" loop provided N D is below a 
certain value. Above this value space charge effects become important and the P — E loops 
will rapidly shrink over a narrow range of Njj until the P r is almost zero. This transition 
happens when the field due to space charge becomes comparable to the depolarising field 
i.e. 



N. 



D 



2P r 
qed 



(15) 



giving the boundary which separates the loops with almost "ideal" from those with negligible 
P r (regions I and III respectively). We should note, that all boundaries in figure^] are not 
sharp, with a continuous but rapid transition between "ideal" and zero P r and E c . 

As depletion widths increase with applied bias a partially depleted sample may become 
fully depleted before the coercive field is reached. In such cases (region II) the hysteresis loop 
will still be similar to that of a sample which is fully depleted at zero bias. In the unshaded 
regions (IV and V) samples are expected to be partially depleted for most biases of interest. 
Within region IV the P — E loops should become slimmer with increasing Njj. Across the 
two blue lines (dashed and solid) the coercive fields drop to 5% and < 1% respectively of the 
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"ideal" value such that in region V the coercive field is expected to be negligible. It may seem 
that regions IV and V are very desirable for applications due to the reduction of E c without 
loss of P r , however in practice increasing Njj also increases the free carrier concentrations 
inside the film and thus the dielectric loss which leads to tilting of the loops as shown by 
Morozovska, as well as to higher leakage currents which bring their own detrimental effects. 
Therefore, experimentally the loops observed in region IV might be either leaky or slim and 
tilted. At this point we should note that our assumption of neglecting the field in the space 
charge free regions of the partially depleted samples becomes invalid when the thickness 
becomes large enough that the resistivity of the bulk becomes comparable to that of the 
depletion layers thus the arguments above will become invalid for thick films beyond the 
dashed curve in figure The position and shape of this curve is difficult to calculate 
(hence the question mark) and in fact is probably very sample specific. For thicker films 
the uniform doping approximation is also not valid as the dopant concentration may vary 
by several orders of magnitude from interface to bulk. 
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FIG. 11: Predicted thickness dependence of coercive field and voltage for Njy = 10 19 cm -3 . The 
critical thickness at which the transition from full to partial depletion occurs will depend on the 
dopant concentration. 

The boundaries in figure ITU1 are only approximate, depending on sample-specific parame- 
ters like P r and e, however, they should give a useful consistency check when interpreting 
experimental data. For example it makes no sense to interpret current-voltage data in terms 
of a partially depleted Schottky model if the hystersis loops are well fitted by the cubic 
equation indicating full depletion; alternatively strongly distorted P — E loops indicate that 
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the field distribution is probably not uniform inside the film. In their recent work [4], Pin- 
tilie et al. have extracted a space charge concentration of 1.8 x 10 19 cm~ 3 for a 250nm PZT 
film. In figure Uni this falls into the partially depleted region IV, consistent with the authors' 
interpretation of the leakage mechanism. 

It is also interesting to look at the thickness dependence of the coercive field for a par- 
ticular dopant concentration. Well within region I the coercive field should be thickness 
independent as shown in figure [TTJ It is then expected to decrease as the boundary is 
crossed and the partially depleted region IV is entered. Within this region polarization 
revesal is initiated by the strong fields within the depletion widths which are thickness in- 
dendent and thus it is the coercive voltage and not field which is expected to be constant. 
Examples of coercive field saturation for low d include the Langmuir-Blodgett polymer film 
data of Ducharme j^J (the interpretation of which, however, has been scrutinised by many 
researchers) and those of Jo et al. on ultrathin BaTi03, which show a constant coercive 
field down to a thickness of 5nm. Thickness independent coercive voltages were reported by 
Trithaveesak et al. [22] for 100-400nm epitaxial BaTiC>3 films and by Kawashima and Ishi- 
wara j^] for 40-65nm SrBi2Ta20g (SBT). Recently Pintilie and Alexe also observed coercive 
voltages independent of thickness (to within 10%) in their epitaxial SrRuOa/PZT/SrRuOs 
structures with d in the range 50-270nm |24J. For thicker films, when the bulk resistivity 
Pt becomes important, switching still requires a constant coercive voltage V c to be dropped 
across the depletion layer, but now the measured coercive field will be given by 



E P . = Vr. 



I _|_ L\ 
d Rd 



where Rd is the depletion layer resistance. This behaviour has been frequently observed for 
thicker films e.g. by Sidorkin j^] in lead titanate and by Merz in 1956 in BaTiC>3 Merz 
argued that the 1/d behaviour was caused by a voltage drop across a surface layer which we 
can identify with a space charge depletion layer as was suggested even earlier by Kanzig [27] . 

Finally we should mention that the highly non-linear field distributions expected for 
higher dopant concetrations have certain implications for conduction in ferroelectrics. In 
the case of fully depleted films the appearance of non-linearity in E(z) coincides with the 
disappearance of the measurable polarization P and thus one is not likely to observe any 
strange conduction bahaviour if these films exhibit ferroelectric P — E loops. However, 
this is not the case for partially depleted films. For capacitors with non-zero Vu the two 
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depletion widths will not be equal even at zero applied bias, as the field will be parallel to 
the polarization in one of them and antiparallel in the other. The depletion region with 
E antiparallel to P will be wider and the maximum field inside it will be lower. As the 
external field is applied to increase the voltage drop across this depletion region (i.e. the 
applied field is trying to switch the film) the magnitude of the field near the interface may 
actually decrease resulting in possible negative differential resistivity if the usual Schottky 
thermionic emission mechanism is responsible for the current flow. Upon exceeding the 
coercive field and switching of the sample, the magnitude of the field at the interface will 
increase discontinuously resulting in a large increase in current. The large difference between 
switching and non-switching currents has been observed previously by Bloom et al 
whereas the negative resistivity is only expected for a certain quite narrow range of and 
thus may not be so common. The large sudden change in resistivity near the coercive field 
can be exploited in resistive switching devices. 



IV. DISCUSSION 



Having presented many interesting results we should now review the assumptions of our 
model and discuss its limitations. Let us start with the gradient term in equation 
which we have chosen to neglect in this work. To our knowledge, the coefficient k has 
not been measured directly. Perhaps the only value so far was obtained by fitting the 
thickness dependence of the Curie temperature to a thermodynamic model, giving k = 
4.5 x 10~ 10 m 3 F~ 1 for lead titanate We have used this value to calculate the relative 
contribution of the gradient term to the free energy and found that even for the highest Np 
considered in our calculations the gradient term amounted to only about 1% of the total 
energy and thus we were justified to ignore it. However, we also ignored surface effects which 
may be quite considerable for very thin films and therefore our description may not be valid 
for very low d. Surface effects have been previously considered by several researchers but 
estimation of their significance is still very difficult due to lack of experimental values for 
parameters such and A. 

Another assumption of our model was that the film is in a monodomain state. It was 
argued by Brennan |8[ that the more stable state may be achieved by relaxation of the 
polarization to form two domains with polarization vectors pointing head to head. However, 
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it was also argued that such a relaxation would require enough charge to be injected to 
compensate the highly charged domain boundary. We have not investigated the possibility 
of domain formation and this should be the subject of future work. 

Finally, we should say a few words about coercive fields. One major drawback of using 
Landau theory for ferroelectrics is that it assumes homogeneous switching and predicts co- 
ercive fields orders of magnitude higher than observed experimentally. Moreover, nucleation 
is known to occur near the interface. The dielectric constant and polarization, however, are 
well described by the thermodynamic approach. We show that the experimentally observed 
low coercive fields can be consistent with Landau theory predictions provided space charge, 
which is always present in reality, is accounted for. 

Nucleation is often observed at specific sites near the interface and many models have 
been proposed to try and undestand how switching occurs. One role of defects or nucleation 
sites was discussed by Ishibashi. In his model, nuclei with reversed polarization are always 
present. The energy barrier for unswitched dipoles is reduced by their interaction with the 
reversed dipoles leading to a reduced coercive field, now dictated by the energy required 
to move a domain wall. Another possibility is the enhancement of local fields by charged 
defects. If this local field exceeds the coercive field a nucleation site will be created. This is 
precisely the effect of space charge. Even for uniformly doped films the field at the interface 
can be orders of magnitude larger than that in the bulk or the measured E = V/d, thus as far 
as nucleation is concerned, space charge acts as an extended defect at the interface. As was 
shown, even moderate dopant concentrations of about 10 19 cm~ 3 can lower the measured 
coercive field by orders of magnitude compared to the intrinsic LGD value. Local space 
charge fluctuations may be responsible for the observed preferential nucleation at specific 
sites rather than homogeneously across the whole interface. 

Whatever the nature of nucleation, domain wall motion or in fact the whole switching 
process, it is clear that the coercive field will, to a large extent, be determined by the local 
field near the site where switching begins and this in turn is obviously heavily influenced by 
space charge. Thus the trends for coercive fields discussed above should be useful at least 
qualitatively. 
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V. CONCLUSIONS 



We have used the phenomenological Landau-Ginzburg-Devonshire theory to investigate 
numerically the effects of inhomogeneous field distributions due to stationary space charges 
on the ferroelectric properties of ferroelectric films. Cases where the samples are fully and 
partially depleted were considered. For fully depleted films both, polarization and coercive 
field values were found to decrease with increasing impurity concentration, as previously 
found by other authors. The field profile in samples with high No was found to have a cubic- 
like spatial distribution, unlike the linear profile expected for semiconductors; this limits the 
range of validity of usual semiconductor models for describing ferroelectrics. Overall the 
effect of space charge was to rescale the LGD coefficients, with the general shape of the 
P — E loops still described by the usual cubic function. 

In the case of partially depleted films, the remanent polarization depends strongly on the 
value of the built-in bias. The P — E loops were found to be more square than the ideal 
space-charge-free loops or the fully depleted loops. A method was proposed which could 
help identify fully depleted films from temperature dependent capacitance and hysteresis 
measurements alone; this may be of particular interest for commercial device characterisa- 
tion, where the physical parameters such as the ferroelectric thickness or composition are 
fixed. 

The unusual field distributions expected for samples with high space charge concentra- 
tions may give rise to interesting current- voltage characteristics such as negative differential 
resistivity. 
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